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Abstract 

Starting from the Hamiltonian operator of the noncompensated two-sublattice 
model of a small antiferromagnetic particle, we derive the effective Lagrangian of 
a biaxial antiferromagnetic particle in an external magnetic field with the help of 
spin-coherent-state path integrals. Two unequal level-shifts induced by tunneling 
through two types of barriers are obtained using the instanton method. The energy 
spectrum is found from Bloch theory regarding the periodic potential as a super- 
lattice. The external magnetic field indeed removes Kramers' degeneracy, however 
a new quenching of the energy splitting depending on the applied magnetic field 
is observed for both integer and half-integer spins due to the quantum interference 
between transitions through two types of barriers. 
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1 Introduction 

The magnetization vector in solids is traditionally viewed as a classical variable. In recent 
years, theoretical and experimental works have demonstrated, however, that the vector 



can tunnel quantum mechanically out of metastable magnetic states or resonate between 
two degenerate ground states[l-13] known as macroscopic quantum phenomena(MQP) 
which are distinguished into the macroscopic quantum tunneling (MQT) and the macro- 
scopic quantum coherence (MQC) respectively. Quantum tunneling of the magnetization 
vector in small single-domain ferromagnetic (FM) particles [1,2], quantum nucleation of 
FM bubbles [3], and quantum depinning of a domain wall in bulk ferromagnets[4] are 
typical examples of the macroscopic quantum phenomena. Similar effect also exists in 
small single-domain antiferromagnetic (AFM) particles in which the Neel vector plays a 
role of macroscopic variable and can tunnel between orientations of lowest energy[5,6]. 
Since the tunneling rate in the AFM particles is much higher than that in the FM par- 
ticles [7], the AFM particles are expected to be a better candidate for the observation 
of MQP than the FM particles. Another interesting phenomenon relating to tunneling 
in magnetization is that for spin systems with discrete rotation symmetry of two folds, 
the tunneling rate is completely suppressed for half-integer total spin known as Kramers' 
degeneracy[8] . Such an effect is called topological quenching in literature [9] and has been 
studied extensively [8- 14] . 

In literature the AFM particle is usually described by the Neel vector of two coUinear 
sublattices whose magnetizations are coupled by strong exchange interaction. External 
magnetic field does not play a role since the net magnetic moment vanishes for idealized 
sublattices. The quantum and classical transitions of the Neel vector in antiferromagnets 
has been well studied[15] in terms of the idealized sublattice model. The temperature de- 
pendence of quantum tunneling was also given for the same model [16] and the theoretical 
result agrees with the experimental observation [17]. A biaxial AFM particle with a small 
non-compensation of sublattices in the absence of external magnetic field was studied in 
Ref.[18] where it was shown that the noncompensated magnetic moment leads to a mod- 
ification of oscillation frequency around the equilibrium orientations of the Neel vector. 
In the present paper we demonstrate that the small noncompensated magnetic moment 
obtains extra energy in magnetic field which changes the original equilibrium orientations 
of the Neel vector and results in interesting tunneling effects in AFM particles. With the 
help of spin-coherent-state path integrals we convert the spin system into a pseudoparticle 



moving in a effective potential —1^(0) with a periodically recurring asymmetric twin bar- 
riers which lead to two kinds of instantons. The total effect of tunnehng gives rise to the 
level splitting which is determined with Bloch theory regarding the periodic potential as 
a superlattice. We show explicitly that the external magnetic field removes the rotation 
symmetry of two-fold and , therefore, the Kramers' degeneracy [10]. The level splitting is 
not quenched any longer for half-integer spin. However a remarkable observation is that 
quantum interference between transitions through two-type of potential barriers results 
in an oscillation of level splitting with the external field. The sphtting could be entirely 
suppressed at the certain value of magnetic field due to the disconstructive interference. 

2 The effective Lagrangian of a biaxial AFM particle 
with a small non-compensated magnetic moment 
and the equilibrium orientations of the Neel vector 

Consider a biaxial AFM particle having two coUinear FM sublattices with a small non- 
compensation.We assume that the particle possesses a x easy axis and x-y easy plane, 
and the magnetic field H is apphed along the y direction. Regarding each sublattice as a 
FM particle the Hamiltonian operator of the AFM particle has the form 

H = E {k^sf + k\\sy^ - ^Hsi) + js, ■ s„ (1) 

a=l,2 

where k±, k\\ > are the anisotropy constants, J is the exchange constant, 7 is the 
gyromagnetic ratio, and §1 , §2 denote the spin operators in two sublattices with the 
commutation relation [5"^, 5"^] — iheijkSaisS^ — x,y,z ; a, /3 — 1,2). Here we 

emphasize that because of the non-compensation of the two collinear FM sublattices the 
interaction terms with magnetic field, i.e. the third term in the summation of Eq.(l), do 
not vanish and result in the equilibrium-orientation change of the Neel vector. We begin 
with the evaluation of the matrix element of the evolution operator in spin- coherent-state 
representation by means of the coherent state path integrals 



(Ny|e-2^H^/'^|N,) = J 



M-l 



n M^k) 

.k=M 



M 



n(Nfe|e--«/'^|Nfe_i) 



.k=l 



(2) 



Here we define |N) = |ni)|n2) , |Nm) = |N/) = |nij)|n2,/) , |No) = \Ni) = |ni,i)|n2,i) 
and tf — ti — 2T , e — 2T/M, respectively. The spin coherent state is defined as 



■\Sa,Sa), (a = 1,2), 



(3) 



where n,^ = (sin 9^ cos <;/)„, sin 9^ sin cos 9^) is the unit vector, = sin (paS^ — cos (paS^ 
and \Sa,Sa) is the reference spin eigenstate. The measure is defined by 
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d/i(Nfc) = Yl — ^- — dn«,fe , dnc,,fe = sin6'a,fcd6'„,fcd0a,fe- 
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As M ^ oo , e ^ 0, exp{—iell/h) ^ 1 — ieR/h, with 



H = E 

a=l,2 

J 
+ - 
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'{S^St + StS^)+2SlSl], (5) 

where •S'^ = 5*^ + ^5*^ , = S^ — iS'a- Making use of the following approximation 

(Nfe|H|Nfe_i) ^ (Nfe|H|Nfe)(N,|N,_i) (6) 



with 



(Nfc|Nfc_i) = Yl (na,fc|na,fc-i), (7) 

a=l,2 

(Nfc|H|Nfe) = E ['^a(^± cos^ ^oc,k + k\\ sin^ ^«,it sin^ (f)a,k) - iHS^ sin ^c«,fe sin 0c.,it] 



a=l,2 



+JS'iS'2[sin 6^1,^ sin 6^2,^ cos(0i,a; - 02,it) + cos6'i,fc cos6'2,fc], 
(na,ik|n„,fe-i) = / 1 + n^.fc n^.fc-i \ exp[-i5'c,Aa(nc,,fe, nc,,fe_i, no)] 



(8) 



exp[-i5'a(0„,fe - 0a,fc_i)(l - cos6'a,fc)], (9) 

and AQ,(ncj jfc, nQ,,jk_i, no) being the area of the spherical triangle with vertices[19] at n^, , n^ j 
and no = (0, 0, 1). Under the large 5" limit we obtain 

(N/le-^^^^/'^lNi) =e-*^«=i.2^"^*«'^-'^"'') / H 2^[0a]2^[^a] exp f'cdt). (10) 

a=l 2 \fl Jti J 

Where £ = £o + is the Lagrangian with 

>Co = E 'S'a^aCOS^a — ■^'S'i5'2[sin^i sin^2cos(02 — 0i) + cos^i cos^2], (11) 

o=l,2 

A = - E (^^'^aCos2e„ + A;||>5^sin2^„ sin2 0„ + 7//5„ sin ^„ sin (/)„). (12) 



For our interest of quantum transition between macroscopic states only the low energy 
trajectories with almost antiparallel Si and S2 contribute to the path integral[18]. We 
therefore replace 62 and (^2 hj 62 = tt — 9i — eg and 02 = + 0i + e^, where ee and 
denote small fluctuations. Working out the fluctuation integrations over ee and the 
transition amplitude Eq.(lO) reduces to 

(N/|e-'"'^/^|N,) = e-^'^^^'^f-'t'i) j V[e]V[4)\ exp Q j'j C'dt^ , (13) 



7 27"' 



V{0,cl>). (14) 



Where (6*1,01) has been replaced by (6*, 0) . 1^(6*, 0) = Q(X_l cos^ 6* + sin^ 6'sin^ - 
mi? sin ^ sin 0) , 5*0 = -^i + 5'2 , m — jh{Si — 82)/^, is the volume of the AFM par- 
ticle. = 2k^S'^/Q and K\\ = 2k\\S^/Q arc the trasnvcrse and longitudinal anisotropy 
constants, respectively. Wc have set Si = S2 = S except in the terms containing — 82- 
The parameter x± — 1^ I ^ is introduced according to Ref.[20] for the problem at hand. 

We consider a very strong transverse anisotropy i.e. K±_ ^ In this case Neel 
vector is forced to lie in the x-y plane. Replacing 6 by 7r/2 + 77^ where r]e denotes the 
small fluctuation and carrying out integral over rjQ we obtain 

(N;|e-"^/^|N,) = e-^°('^/-^^) j P[0] exp C^jfdt) , (15) 

where 

A// = (// + /») (16) 

is the effective Lagrangian which is seen to be the Lagrangian of a plane rotor. Where 
If — m^/27^ir_L and = X^ll^ ^'^^ effective FM and AFM moments of inertia, 
respectively [6]. T^(0) = Jl-fr||(sin0 — A)^ (A = H/Hc with a parameter He — 2K\\/m) is 
the effective potential. It is seen that the net magnetic moment of the noncompensated 
sublattices in the applied magnetic field shifts the equilibrium orientations of Neel vector 
for corresponding angles iarcsinA as shown in Fig.l-(b) besides the modification of 
FM moment of inertia If given in Refs.[6,18]. It may be worth while to compare our 
results with that in literature. In the absence of the magnetic field (namely A = 0) the 
two degenerate equilibrium orientations return to the positive and negative x-axis (see 
Fig.l-(a)) respectively in agreement with the equilibrium phases of the AFM particle with 



noncompensated sublattices[18]. The small oscillation frequency of Neel vector around 
its equilibrium orientations which serves as a characteristic parameter for the flip of Neel 
vector of the AFM particle is seen to be 



u}[H, m 

For A = it reduces exactly to 



u{H = 0,m) 



h + h 



2K\\ 



(17) 
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as given in Ref.[18]. If we consider idealized sublattices that m = the frequency Eq.(17) 
goes back to the well known value 



uj{H — m — Q) — ^ 



2Kii 
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2KnJ 
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(19) 



The effective potential of the plane rotor is plotted in Fig. 2. The minima of the potential 
correspond to the equilibrium orientations of the Neel vector. The energy of net magnetic 
moment in the applied magnetic field lowers the barrier height in the direction of magnetic 
field while increases the barrier height in the opposite direction. We are interested in the 
quantum tunneling of the effective plane rotor through the barriers. 



3 Two types of instantons and level shifts 

In order to obtain the tunneling rate we evaluate the Euclidean path integrls in Eq.(15) 
with the Wick rotation t = ir . The Euclidean Lagrangian for the pseudoparticle 
moving in the classical forbidden region, namely, in the barrier is seen to be = 
{If + Ig)^ (^^^ + V{(l)). The equation of motion of rotor at finite energy is 

The periodic potential V{(j)) — V{(f) + 2n7r) has an asymmetric twin barrier (Fig.2) .The 
two-fold rotation symmetry[10], namely, y(0-|-7r) = V{(f)) in the presence of magnetic field 
is removed . When the energy is higher than the ground sta8te tunneling is dominated by 
periodic instantons[21,22]. Thus there are two different periodic instantons corresponding 



to two types of barriers. With the periodic boundary condition, two periodic instanton 
solutions of Eq.(17) are found to be 

= - + 2 arctan[Aisn(gT, k)] , (21) 
2 

(f)^^^ = 2 arctan[A2sn(gr, k)] , (22) 



where sn{qT,k) is eUiptic function with modulus k and period AK{k). K{k) denotes the 
complete elliptic integral of the first kind, with 



k 



1/2 



The parameters q , Ai and A2 are defined by 
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nl/2 



[1 - (-1)^A)2 -772 



, (i = l,2) 



.2(// + /a), 

The trajectories of instantons 0^^^ and 0^^^ are shown in Fig.2. At initial time r^, instantons 
(f)'^^^ , 0^2^ start from the potential well at 0, = arcsin A and reach the neighbouring well 
at (f)f = TT — arcsin A at final time t/ along the anticlockwise (through small barrier) 
and clockwise (through large barrier) paths respectively. In other words the Neel vector 
tunnels through a large barrier (or small barrier) between two angular positions with the 
lowest energy. 



We assume that |m, 02n) and |m, 02n+i) denote the eigenstates of the harmonic os- 
cillator approximated Hamiltonian in the potential wells at 02n = 2n7r + arcsin A and 
4>2n+i = (2?^ + l)'7r — arcsin A, respectively, where m is the index of low- lying levels. The 
amplitudes tunneling through two different barriers are given by [21] 
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(23) 



(24) 



where A£^)(A£^)) is the level shift induced by tunneling through the small (large) barrier 
alone. AgW is actually the overlap integral defined in the following Eqs.(42),(43). Where 
P = Tf — Ti. With the help of the path integral, the matrix element in Eq.(23) and Eq.(24) 
can be rewritten as 

= / '/•/)^m(</'!'\ <l>i)}C^'\h, rf, 0., rOd0/d0, , (25) 



where 



= / C(0?\</'/)V'm(0^'\0O^^'H0/,T/; 0^,T,)d0/d(/), , (26) 



c(0 
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1,2). 



(27) 



(28) 



Here S^^' denotes the Euchdean action, and K,^'^\(t)f,Tf\ (pi,Ti) is Feynman propagator 
through two kinds of barriers. Substituting the periodic instanton solutions Eq.(21) and 
Eq.(22) into Eq.(28) the Euchdean action along the classical trajectory is found to be 



X^^E{k) + (fc^ - Xf)Kik) + (Xf - fc2)n(fc. A,) 



(29) 
(30) 



where E{k) is the complete elliptic integral of the second kind, and n(A;, Aj) is the complete 
elliptic integral of the third kind with the parameter Aj. The level shift AsW can be 
determined by completing the integrals of (25) and (26) and comparing the result with 
Eq.(23) and Eq.(24). Using the method in Ref.[21], the transition amplitude is obtained 



as 



exp 



sinh 
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2Kii 
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1,2), (31) 



where k' = (1 — k^y/'^ and a = [(1 + 77)^ — A^]-^/^. Comparing this expression with Eq.(23) 



and Eq. (24) we find that two types of level shifts are 
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exp 



1,2), 



(32) 



4 The magnetic field dependence of tunneling rates 

It is easy to see that the difference between the heights of larger and small barriers 
increases with the external magnetic field . On the other hand, the effective frequency 



of oscillator near the bottom of potential well given in Eq.(17) i.e. u — ujq{1 — A^)^/^ 
decreases with the increasing field, where uq — {j-j^Y^'^ is the frequency in the absence 
of the magnetic field[18]. For the low energy case that the energy E is far below the 
barrier height , i.e. <^ 1 — A, /c ^ 1, /c' — 0, E{k:), K{k) can be expanded as power 
series of k' . The complete elliptic integral of the third kind is expressed as 

nl/2 



n(A;,A, 



+ A? 



{E{k)F{a\k')+K{k)[E{a\k')-F{a\k')]}, 



(l + Af)(F + Af)J 

where = arcsin(A^/(/c^ + A^))^/^. F(q;(*\A;'), E{a^''\k') are the incomplete elliptic 
integrals of the first and second kinds, respectively. Thus ^{k, Aj) can be also expanded 
as the power series of k' . Then, the power series expression of W^'^^ reads 



jy(i) ^ '^ML J(i _ a2)1/2 _ Aarecos a - \(1 - 
ojq I 16 

^(2) ^ 4^ r ^ _ ^2)1/2 ^ ^ arccos(-A) -^{l- A'f'k" [in ^ + 11) 
iOf] L Id l k J J 



(33) 



(34) 



In the low energy case, k' = 4?7/(l — A^) <C l,we may take the oscillator approximated 
energy-quantization i.e. E E^ — {'m+ l/2)tiuj. Taking note of limits K{k ^ 0) ^ 7r/2 
and (7(77 ^ 0) ^ (1 - A2)-V2 find 



AsW = F^(l - A2)5/4+3W2 exp{-S[(l - A')V2 _ Aarecos A]}, 
A£(^) = Fm{l - A2)5/4+3m/2 exp{-B[(l - A^) V2 + A arccos(-A)]}. 



(35) 
(36) 



Where 



huJn 



TT 



1/2 



Eq.(35) and Eq.(36) give rise to the field dependence of the level shift for low- lying ex- 
cited states. There is an obvious difference between the level shifts induced by tunnehng 
through two kinds of barriers. For a given excited state, Ae^^^} as a function of the external 
magnetic field is plotted in Fig. 3, with k\\ = lO^erg/cm^ , K± = lO^erg/cm^ , x± = 10^ ? 
the excess of spin S1 — S2 — 10 and the AFM particle radius r = l.hnm. It is clearly shown 
that the tunneling rate through a small barrier increases rapidly with the external mag- 
netic field because the field reduces both the height and width of barrier. The situation is 
just opposite for the tunneling through the larger barrier . In addition. Fig. 3 also shows 



that the tunnehng rate increases with the energy levels in the low-lying excited states. 
When A = H/ He attends to 1 there is no tunnehng at all since the small barrier shrinks 
to zero and only one easy direction remains. In the absence of applied magnetic field we 
have (A = 0)A£:^-' = Ae^-* and for the ground state tunneling,namely E = rj = 0,the 
level shift Asq reduces to exactly the result in Ref.[18]. 



5 Level splitting and quantum interference effect 

Ae^^ is only the level shift induced by tunneling through a single barrier (smaller or 
larger). The periodic potential V{(j)) = V{4> + 2n7r) can be regarded as a one-dimensional 
superlattice consisting of two sublattices. The general translation symmetry results in 
the energy band structure, and the energy spectrum could be determined with the Bloch 
theory. Let \m, 02n) be the eigenstates of the zero order Hamiltonian Hq^^ in the potential 
well which lies at 02n = 2n7r + arcsin A. 02n+i) denote the eigenstates of the zero 

(2) (2) 

order Hamiltonian Hq in the well at 4>2n+i = (2^ + l)7r — arcsin A. Thus 

hJ^V^S+i) = £m|m,0S+i). (38) 
Bloch state with 2n periodic boundary condition is written as 

m = E (e^(^+"°)^- \m, 0S) + e^(«+^")^S+. ^^)^ , (39) 

where ^ is Bloch wave vector , 1) is seen to be the topological phase from 

Eq.(15). Substituting Eq.(39) into the following stationary Schrodinger equation 

HIV') = E\i;) , (40) 

and taking into account only the nearest neighbours yield the energy spectrum as 

E ^ Sm- Ae^^) cos[(C + ^o) (tt - 2 arcsin A)] 

- A£(^) cos[(e + So){Tr + 2 arcsin A)], (41) 



where the level shift A^W is actually the overlap integral defined by 

A£(^) = -/<(0-0S)H«„.(0-0g+i)d0, (42) 
A£(^) = - 1 <(0 - 0S)Hi.^(0 - 0g_Od0 . (43) 

The Bloch wave vector ^ can be assumed to take either of the two values and 1 in the 
first Brillouin zone[21,23]. Thus the level sphtting is seen to be 

Aem = I As^^) cos[2(l + So) arccos A] + (-1)^^° As^^ cos[2(l + So) arccos A] 
-Ae(^) cos[2^o arccos A] - {-1)^'^'> As^^^ cosp^o arccos A] | 
= £;^|sin[(2S'o + l) arccos A] I , (44) 

where 

£■+ = i?^cosh(B7rA/2) , (^o = integral) , (45) 
E- = i?^sinh(S7rA/2) , (^o = half - integer) , (46) 

= 4F^(1 - A2)^/^+3W2 

exp{-B[(l - A^)^/^ + Aarcsin A]}. (47) 

When H = 0, i.e. arcsin A = 0, energy spectrum in Eq.(41) reduces to the re- 
sult in Ref.[14], and for 5*0 = half-integer the MQC is quenched in agreement with 
Kramers' theorem which can be well understood as two-fold discrete rotation-symmetry 
of Hamiltonian[10] . The apphed magnetic field breaks the rotation symmetry and the 
Kramers' degeneracy is removed. Level splitting Eq.(44) shows the quantum interfer- 
ence effect depending on applied magnetic field. The level splitting increases with the 
magnetic field and whenever the magnetic field reaches some specific values which satisfy 
H/Hc — cos[Z7r/(25'o -I- 1)] {I is a integer), we have As^ = no matter 5*0 is integer or 
half-integer. The quenching is similar to the case of the FM particle in Refs.[9,24] and is 
the result of quantum interference between transitions through the two kinds of barriers. 
Fig. 4 shows the oscillation of the level splitting with the field for the ground state. 



6 Conclusion 



We present a full study of quantum tunneling effect for AFM particle with a small non- 
compensation of sublattices in an external magnetic field. The level splitting which is 
obtained only for ground state in literature has been extended to low-lying excited states 
with the help of periodic instanton method. The quantum interference effect, particularly, 
the entire suppression of tunneling may be of significance for practical application. 
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Figure captions 



Fig.l (a) The two degenerate equilibrium orientations of Neel vector with noncom- 
pensated sublattices in an apphed magnetic field. 0+ = arcsinA. (b) The equilibrium 
orientatons of Neel vector in the absence of applied magnetic field (A = 0). 

Fig. 2 The periodic potential with asymmetric twin-barrier and the instanton trajec- 
tories. 

Fig. 3 The level shift AsW as a function of H/Hc- Sohd line for AeW and dotted fine 
for Ae(^\ 

Fig. 4 The level sphtting at ground state as a function of H/ He- 
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